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Abstract 

This paper deals with the reconstruction of the direction of a gravitational wave source using the 

(N 

' detection made by a network of interferometric detectors, mainly the LIGO and Virgo detectors. 

00 : 

We suppose that an event has been seen in coincidence using a filter applied on the three detector 
, data streams. Using the arrival time (and its associated error) of the gravitational signal in each 

\ detector, the direction of the source in the sky is computed using a minimization technique. 

O^. Por reasonably large signals (SNR>4.5 in all detectors), the mean angular error between the real 

. location and the reconstructed one is about 1°. We also investigate the effect of the network 

bp. 

geometry assuming the same angular response for all interferometric detectors. It appears that 
^ ' the reconstruction quality is not uniform over the sky and is degraded when the source approaches 

the plane defined by the three detectors. Adding at least one other detector to the LIGO- Virgo 
network reduces the blind regions and in the case of 6 detectors, a precision less than 1° on the 
source direction can be reached for 99% of the sky. 
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I. INTRODUCTION 



The LIGO and Virgo gravitational wave interferometric detectors are approaching their 
design sensitivity Q],!^] and in the near future, coincidences between the three detectors 
(LIGO-Hanford, LIGO-Livingston and Virgo) will be possible. In order to reconstruct the 
direction of the astrophysical sources in the sky, it is well known that a minimum of three 
detectors is mandatory even if an ambiguity remains between two positions symmetric with 
respect to the plane defined by the 3 detectors. The source direction is also provided by the 
cche.en. .a.c.e. fo. S, g o. coal^c.g b.an. Q, 3. B w,e. o„e of .he 

outputs of the detection algorithm is an estimation of the source direction. 

In this paper, we propose a method for estimating the source position using only the 
arrival time of the gravitational signal in each detector. The event detection is supposed to 
have been previously done by dedicated algorithms ([10-23] for bursts, [24-29] for coalescing 
binaries) and is not within the scope of this article. 

The direction reconstruction is based on a minimization as described in Section II. 
This technique can be easily extended to any set of detectors. Moreover, the method can 
be applied to several types of sources (burst, coalescence of binary objects ...) as soon as an 
arrival time can be defined for the event. 

Section II also deals with the simulation procedure which will be used in the following 
sections to evaluate the reconstruction quality in several configurations. Sections III and IV 
describe the performances of the LIGO- Virgo network first neglecting (III), then including 
(IV) the angular response of the detectors. In section V, we consider the addition of other 
gravitational wave detectors (supposing a similar sensitivity) and investigate their impact 
on the reconstruction. In real conditions, systematic errors on arrival time are likely to exist 
and their impact on the reconstruction is tackled in Section VI. 



II. MODELING THE RECONSTRUCTION OF THE SOURCE DIRECTION 
A. Arrival time of GW signals in interferometers 

Within a network of n interferometers, we suppose that each detector Di measures the 
arrival time ti of the gravitational wave. Of course, the definition of the arrival time depends 
on the source type and is a matter of convention, for example: peak value in the case of a 
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supernova signal, end of the coalescence for binary events. In the following, it is assumed 
that all interferometers use consistent conventions. 

The error on the arrival time, (jj, depends on the estimator used and on the strength of 
the signal in the detector Di, strength (for a given distance and a given signal type) which is 
related to the antenna pattern functions (see |5i|] and references therein) at time tj. At that 
time, the antenna pattern depends on the longitude and the latitude of the detector location, 
as well as its orientation, the angle between the interferometer arms, the sky coordinates a 
(Right Ascension) and 5 (Declination) of the source, and the wave polarization angle ip. 

The timing uncertainty cTj can be parametrized by jlfl| : 

{SNR{)<^ ^ > 

where SNRi is the measured SNR in detector Dj, ctq and ^ are constants depending on 

the detection algorithm and the signal shape. For example, a burst search with a 1-ms 

Gaussian correlator leads to ctq ~ 1.4 ms and C = 1- Typically, for an SNR equal to 10, the 

error on the arrival time is a few tenth of milliseconds and weakly depends on Q for SNR 

values between 4 and 10. 



B. Definition 

The n measured arrival times tj and their associated errors cTj are the input for the 
reconstruction of the source direction in the sky, direction defined by a and 5. 

In the 3-detector configuration, the angles {9 and ) of the source in the detector coor- 
dinate system (see Ref.Q| and 3] for exact definitions) are given by 0: 
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or a least-square function involvinsr 

n 

^. This minimization also implies 



where Di is placed at (0,0,0), D2 at (a, 0,0) and D3 at (61,62,0). 

When performing a coherent analysis of the GW detector streams the position of the 
source in the sky is part of the output parameters, corresponding to the stream combination 
which maximizes the SNR. However, for a burst searchjjt is known that thousands of possible 
positions have to be tested to obtain the solution 
the integration of detector streams has to be minimized 
the test of hundreds of initial conditions in order to reach the right minimum. 

Concerning coalescing binaries, it implies the definition of a five-parameter bank of filters 
including the chirp mass, the three Euler angles and the inclination angle [3] of the orbital 
plane or a three-parameter bank of thousands filters for the two source angles and the chirp 
time ^. 

In all coherent techniques, the extraction of the source direction is an heavy process 
imbedded in the detection procedure. 

In this paper, we propose a simpler approach where a and 5 are found through a least- 
square minimization using separately triggered events obtained by a coincidence search. We 
suppose that the detection is already performed applying suitable algorithms (matched filter 
for coalescing binaries, robust filters for bursts). The is defined by: 

, ^ (t.-(to + Af-^^(c.,^)))^ 

where to is the arrival time of the gravitational wave at the center of the Earth and 
Af (a, 5) is the delay between the center of the Earth and the i*'' detector which only 
depends on a and b. 

The first advantage of this definition is that it deals with absolute times recorded by each 
detector rather than time differences where one detector has to be singled out. Otherwise, 
the best choice for the reference detector is not obvious: the detector with the lower error 
on the arrival time, the detector which gives the larger time delays or the detector leading 
the best relative errors on timing differences ? This definition leads to uncorrelated errors 
on fitted measurements. 

The second advantage is that the network can be extended to any number of detectors 
and the addition of other detectors is straightforward. Obviously, the method requires that 
the event is seen by all detectors. 
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The least-square minimization provides the estimation of to,a,S and the covariance ma- 
trix of the fitted parameters. When the number of detectors is greater than 3, the value at 
the minimum can also be used as a discriminating variable, as the system is overconstrained. 



C. Simulation procedure for a 3-detector network 

A list of coincident events are defined by the three arrival times and their associated 
errors. No detection procedure is performed in these simulations as stated before. 

The simulation proceeds in two steps in order to study two coupled effects: antenna- 
patterns and location with respect to the 3-detector plane. The first step is a simplified 
approach: the antenna-pattern functions are ignored and the same error ai is assumed for 
arrival times. The second step is more realistic: we assume the same sensitivity for each 
detector and the signal strength is adjusted to have the mean (over the three detectors) SNR 
equal to 10. However, this implies that sometimes, due to the antenna-pattern, the signal is 
seen in a given detector with an SNR lower than 4.5, which remains an acceptable threshold 
for a real detection. The same threshold equal to 4.5 will be used later to see the effect of a 
reasonable detection scheme on the angular reconstruction error. We evaluate the errors ai 
for each arrival time using Equation ^ with C = 1- 

For a given simulation, the coordinates (a*™^, 5*''"^) of the source are chosen and the true 
arrival times tf"^'^ on each detector are computed taken tg''"^ equal to (it is obvious that 
the timing origin can always been chosen such as tg''''*^ = 0) 

The measured arrival times t'^^"-^'^'^^'^ are drawn according to a Gaussian distribution 
centered on t*™^ and of width a^. The simulated values t'^^"'^^^'^'^ are then used as inputs for 
the least-square minimization. As the minimization is an iterative procedure, some initial 
values for the parameters have to be given, to is initialized by the average of t™'^'^*"''^'^. For the 
angles, it appears that the initial values for the angles have no influence on the minimization 
convergence and a random direction is adequate. In the case of three interferometers, it is 
well known that there is a twofold ambiguity for the direction in the sky which can lead to 
the same arrival times in the detectors. These two solutions are symmetric with respect to 
the 3-detector plane. In order to resolve the ambiguity, a fourth detector is needed. For the 
evaluation of the reconstruction accuracy, only the solution closest to the source is retained. 
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FIG. 1: Reconstruction accuracy for a source at the Galactic Genter at a fixed day time for LIGO- 
Virgo network. 

The angular error is the angular distance on the sphere between the true direction and the recon- 
structed one. The curves on the first two plots correspond to the best Gaussian fit. 

III. SOURCE DIRECTION RECONSTRUCTION: 3-DETECTOR NETWORK 
GEOMETRY 

In this section, we only deal with the LIGO- Virgo network and the effect of the antenna- 
pattern functions are not included and it is assumed that all detectors measure the arrival 
time with the same precision. As previously said, it allows to decouple the effect of the 
antenna-patterns and of the location with respect to the 3-detector plane. 

First of all, as an example, in order to evaluate the accuracy of the reconstruction, we 
choose a given position in the sky (coordinates of the Galactic Center ace = 266.4°, 6gc = 
—28.98°) and we perform the simulation with cTj = 10^^ s at a fixed time. The results are 
shown on Figure ^ A resolution of about 0.7 degrees can be achieved both on a and 6. The 
angular error is defined as the angular distance on the sphere between the true direction and 
the reconstructed one (it does not depend on the coordinate system and in particular there 
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is no divergence (only due to the coordinate system) when S is equal to 90 degrees). This 
variable will be used in the following steps as the estimator of the reconstruction quality. 
The mean angular error is 0.8 degrees. As shown on Table El the estimated errors (given by 
the covariance matrix) obtained by the minimization are in perfect agreement with these 
resolutions. 



Angle 


RMS of Distribution 


Mean of estimated errors given 
by the covariance matrix 


a 


0.760° 


0.758° 


6 


0.635° 


0.636° 



TABLE I: Reconstruction accuracies on a and 5 and errors given by the covariance matrix. Three 
digits are given in order to show the adequacy between RMS and errors given by the covariance 
matrix. 

Varying cr, in the simulation, it is found that the resolution is proportional to the timing 
resolution for reasonable errors (0"^ < 3 ms. For (Xj ^ 3 ms, the errors become comparable 
to the sky size.) and we obtain: 



OiResolution — 0.7° (5) 

r u s 

^Resolution 0.6 Z 
10"*S 

Mean Angular error = 0.8° — ^ 
^ lO-^s 

For the Galactic Center direction and with cTj = 0.1 ms, we evaluate the reconstruction 
performances over one entire day, still neglecting the antenna-pattern. 

It clearly appears on Figure |21 that the resolution varies, as expected, during the day and 
lies between 0.8° and 4.3° (see Table HTl for details). 

The error maximum which appears around t = 20h (the time origin is arbitrary) cor- 
responds to directions for which the source approaches the detector plane. The angular 
distance between the two possible solutions (the source and its mirror image) gives a good 
estimator of this closeness (see Fig. This effect is even clearer on Figure El where the 
source is located at 6 = 0°. In this case, the source crosses the 3-detector plane twice a day 
{t = 2.4 h.,t = 14.4 h) and the angular distance between the real source and its mirror image 
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FIG. 2: Angular Error for a source at the Galactic Center for one day (t=0 arbitrary chosen) for 
LIGO- Virgo network. 

The top-left plot presents the angular error as a function of time. The bottom-left plot shows 
the angular distance between the source (S) and its mirror image (M) (other direction in the sky 
which gives the same time delays). The top-right plot is the distribution of angular errors. The 
bottom-right gives the fraction of events with an angular error below a given threshold versus this 
threshold, ai is set to 0.1 ms for all detectors. 

is null. This degradation of the angular reconstruction is not related to the minimization 
technique. It is a geometrical property of the network which releases constraints when the 
source belongs to the 3-detector plane. For the 6 = 0° configuration, angular resolutions as 
given in Table ITTl are similar to the ones obtained for the Galactic Center. 

This degradation can be easily understood with only two detectors located at {±d/2, 0, 0) 
and a source in the (x,y) plane defined by its angle 6 with the x axis. In this 2-detector 
case, the arrival time difference ^21 is given by: 

d 

t2i = — cos 9 (6) 
and thus, if toi is measured with an error 6t, it will induce an error 59 = ,|^!^* , on the 
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6n 


Minimal Error (°) 


Maximal Error (°) 


Mean Error (°) 


Median Error (°) 


-28.98 (GC) 


0.8 


4.3 


1.9 


0.95 





1.3 


3.1 


1.8 


1.5 


All values 


0.7 


4.5 


1.6 


1.1 



TABLE II: Angular resolution for Galactic Center, 6 = 0° position and averaged on all possible 
values of 6 taken over the whole day. No antenna-pattern effect is included. 

source direction which diverges when the source belongs to the 2-detector line = or tt). 
Of course, this simple estimation of the angular error supposes that Eq. IHl between ^21 and 
9 can always be inverted, assumption which fails when ^21 becomes greater than ^ due to 
measurement errors. We can wonder how to handle this case. In the 2-detector example, 
we minimize defined by: 

X [d) = — U21 - -cos^ 1 (7) 

which becomes minimal for = acos(t^'^"*"'"'^'^ x |) if x^|<lor6' = when 

the previous condition is not satisfied. We observe similar effect in the 3-detector case when 
the source approaches the 3-detector plane. 

Figure ^ presents the angular errors averaged on all possible value of 6. In this case, the 
angular resolutions are given in Table m ranging 0.7 to 4.5 degrees. Without considering 
antenna-pattern effect, it means that an angular error lower than about one degree can be 
reached for half of the sky and 90 % is below 3.5 degrees. The intersection of the 3-detector 
with the celestial sphere appears as a zone with worst angular resolutions (error ranging 
from 3 degrees to 4.5 degrees) of the source direction. 

IV. INCLUDING THE ANTENNA-PATTERN EFFECT IN A 3-DETECTOR 
NETWORK 

As in previous section, we use the LIGO-Virgo network as benchmark. Now, the strength 
of the signal seen by each detector is computed taking into account the antenna-pattern 
functions. The errors on arrival times are estimated using Eq^with C = 1 and o"o = 1ms. 
In order to remain close to the case described in Section III (all timing errors equal to 
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Time of the day (hrs) Angular Error Threshold (°) 



FIG. 3: Angular error for a source at 5 = 0° for one day (t=0 arbitrary chosen) for LIGO- Virgo 
network. 

The top-left plot presents the angular error as a function of time. The bottom-left plot shows 
the angular distance between the source (S) and its mirror image (M) (other direction in the sky 
which gives the same time delays). The top-right plot is the distribution of angular errors. The 
bottom-right gives the fraction of events with an angular error below a given threshold versus this 
threshold, ai is set to 0.1 ms for all detectors. 

0.1 ms), we impose that the mean SNR over the three detectors is equal to 10. A linear 
polarization has been assumed for the incoming signal. 

Figure El presents the angular error as a function of the day time for a source located at 
6 = 0° and can be compared to Fig. El The two broad peaks (t = 2.4 hours and t = 14.4 
hours), corresponding to the time when the source belongs to the 3-detector plane, are still 
present (barely seen due to the change of scale) but six sharp diverging peaks show up. They 
correspond to blind regions for at least one detector, as expected. 

Of course, in real conditions, a minimal threshold will be applied for the event selection. 
For this purpose, we define a SNR threshold equal to 4.5 on each detector which leads 
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FIG. 4: Reconstruction resolution for the whole sky for LIGO- Virgo network. 
The top plot presents the angular error as a function of time for all possible values of S. The 
bottom-left plot is the distribution of angular errors. The bottom-right gives the fraction of events 
with an angular error below a given threshold versus this threshold. The white line gives the 
trajectory of the Galactic Center during the day. cjj is set to 0.1 ms for all detectors. 

(supposing a Gaussian noise) to a false alarm rate about 10~^ Hz when requiring a triple 
coincidence between Virgo, LIGO-Hanford and LIGO-Livingston. About 79 % of the events 
satisfy this constraint and a median error of 2.0° is reached (See Table ITTTl and Fig. Ej). 

For a source located at the Galactic Center, the angular error shows several spikes related 
to the antenna-pattern effect and the geometry of the network only increases the mean error 
in the region around t = 21 h (Fig. [7j). All peaks correspond to regions which are blind for 
at least one interferometer. Imposing that all SNR are above 4.5 leaves 56.4 % of the events 
and in particular all these blind regions are removed. The obtained angular resolutions are 
given in Table IIIIl They are slightly better than for the 6 = 0° case. 

Figure |H1 presents the angular errors averaged on all possible 6. As in Fig. IH the in- 
tersection between the 3-detector plane and the celestial sphere is visible. The regions of 
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FIG. 5: Angular Error for a source at 5 = 0° as a function of the day time (t=0 arbitrary chosen) 
for LIGO- Virgo network. 

The top plot presents the angular error as a function of time. The middle plot shows the angular 
distance between the source and its mirror image (the other sky direction giving the same time 
delays). The SNR values in each interferometer are plotted on the bottom figure. The effect of the 
antenna-pattern functions is included and the mean SNR is set to 10. 

largest errors (> 15°) corresponds to the blind regions of the various detectors. Without 
any conditions on the SNR in each detector, we obtain a median error of 1.7° (see Table 
IIII|) . The SNR condition is satisfied for 60% of the events and leads to a median error of 
1.3° with 90 % of the events below 3.7°. In the most favorable cases, we can even reach a 
precision of 0.7°. A resolution better than 1° is obtained in 30% of the cases. 

V. ADDITION OF OTHER GRAVITATIONAL WAVE DETECTORS 

First of all, the 2-kilometer long Hanford interferometer has been included in the network. 
In this case, we supposed that the SNR seen by by Hanford 2k is half of the SNR seen by 
Hanford 4k. As expected, adding Hanford 2k does not change significantly the resolutions 
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FIG. 6: Distribution of angular errors for a source at 5 = 0° for LIGO- Virgo network. 
The top-left plot shows the distribution of angular error for all events. The bottom-left one presents 
the fraction of events with an angular error below a given threshold versus this threshold. The 
right plots are the same imposing that the SNR seen by each detector is greater than 4.5. 

obtained in previous sections. 

Then, we add TAMA 0, GEO H and AIGO 3 assuming (for sake of simplicity) the 
same sensitivities as LIGO and Virgo. We still impose a mean SNR equal to 10. Figure 13 
which has to be compared to Figure El shows the angular reconstruction for a source located 
at 5 = 0° as a function of the day time adding only a fourth detector to the LIGO- Virgo 
network. Of course, there is no longer an ambiguity in the possible solutions and all former 
blind regions are attenuated except the one around t = 8.5h for which effective low SNRs 
are obtained with both Hanford and Virgo. The mean error (3° for LIGO- Virgo Table 1111)) is 
lowered to 2.0° adding GEO, to 1.7° adding TAMA and 0.7° adding AIGO. Not surprisingly, 
the larger the network baseline, the better the resolution. 

In order to evaluate the performances of a strongly overconstrained network, we evaluate 
the reconstruction resolution for the full sky using the 6 detectors (still using the same 
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6n 


Event Fraction(%) 


Min. Error (°) 


Mean Error (°) 


Median Error (°) 


-28.98 (GC) 
Any SNR 


100 


0.7 


4.0 


1.8 


-28.98 (GC) 
All SNR > 4.5 


56.4 


0.7 


1.8 


1.25 




Any SNR 


100 


1.2 


3. 


2.2 




All SNR > 4.5 


79.1 


1.2 


2.1 


2. 


All values 
Any SNR 


100 


0.7 


2.7 


1.7 


All values 
All SNR > 4.5 


59.8 


0.7 


1.8 


1.3 



TABLE III: Angular resolution for Galactic Center, 6 = 0° position and averaged on all possible 
values of 6 taken over the whole day. The antenna-pattern effect is included. The column "Event 
Fraction" gives the fraction of events satisfying the SNR selection criterion. 

sensitivity). The results are shown in Figure ITUl similar to Figure |H1 The mean error is 
about 0.4° and 99% of the sky is covered with a resolution below 0.8°. The regions with 
errors between 1 and 3 degrees correspond to directions for which several interferometers 
have very low SNRs (as the region around t = 8.5h in the previous figure). 

VI. EFFECTS OF SYSTEMATIC ERRORS ON ARRIVAL TIMES 

All angular errors quoted previously suppose that arrival time measurements are only 
subject to Gaussian noise. Systematic biases can also be introduced by the analysis and 
their effect can be evaluated. In order to do so, we modify Equation |H1 introducing a timing 
bias for only one detector: 

^Measured ^ ^true ^ QaussianRandom xai + Bias (8) 
As in sections lllll and IIVI we only consider the LIGO- Virgo network. It appears that 
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FIG. 7: Angular Error for a source at the Galactic Genter as a function of the day time (t=0 
arbitrary chosen) for LIGO- Virgo network. 

The top plot presents the angular error as a function of time. The middle plot shows the angular 
distance between the source and its mirror image (the other sky direction giving the same time 
delays). The values of the antenna-pattern functions are plotted on the bottom figure. The effect 
of the antenna-pattern functions is included and the mean SNR is set to 10. 



Bias (ms) 


Bias on a{°) 


Bias on 6{°) 


Angular Bias (°) 


.1 


-0.24 


0.65 


0.68 


.2 


-0.48 


1.26 


1.33 


.3 


-0.72 


1.79 


1.90 


.4 


-0.97 


2.55 


2.69 


.5 


-1.2 


3.17 


3.35 


1. 


-2.6 


6.20 


6.63 



TABLE IV: Effect of a timing bias on angular reconstruction. 

The bias is applied to the arrival time at Livingstone interferometer. The RMS of the statistical 
errors has been set to .1 ms. For this direction, the statistical angular error is about 0.8°. 
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FIG. 8: Reconstruction resolution for the full sky for LIGO- Virgo network. The top graph gives 
the angular error as a function of a and sin((5). For clarity, the values above 20° have been set to 
20°. The bottom plots presents the angular error distribution and the fraction of events below a 
given angular threshold imposing or not the SNR condition (all SNR>4.5). 

the widths of the distribution for reconstructed a and 5 are not modified by the bias but 
the central values are shifted from the true ones. The differences between the reconstructed 
value and the true one are proportional to the bias and are significantly different from zero 
when the bias and the statistical error have the same order of magnitude. Table HVl shows the 
effect of the bias for a given direction (we check that the effect is independant of the source 
location). In this example, the bias has been applied to the Livingstone interferometer. The 
width of statistical errors on arrival time was .1 ms leading to a statistical angular error 
about 0.8°. For the tested configurations, we do not observe significant differences between 
the three interferometers of the network. 
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FIG. 9: Reconstruction resolution for a source at (5 = 0° as a function of the day time (t=0 arbitrary 
chosen) adding a fourth detector to the LIGO- Virgo network. 

The top figure presents the angular error as a function of time adding TAMA, the middle one 
adding GEO and the bottom one adding AIGO. 

VII. CONCLUSION 

We described a method for the reconstruction of the source direction using the timing 
information (arrival time and associated error) delivered by gravitational wave detectors 
such as LIGO and Virgo. The reconstruction is performed using a least-square minimization 
which allows to retrieve the angular position of the source and the arrival time at the center 
of the Earth. The minimization also gives an estimation of errors and correlations on fitted 
variables. For a given position, the angular error is proportional to the timing resolution 
and the systematic errors (if they exist) introduce a significant bias on reconstructed angles 
when they reach the level of the statistical one. 

When the antenna-pattern effect is included and imposing a mean SNR value of 10 in 
the LIGO- Virgo network, a precision of 1.7° can be reached for half of the sky. In order to 
reproduce a realistic case, we apply a threshold on the SNR in each detector (SNR> 4.5 
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FIG. 10: Reconstruction resolution for the full sky with a network of 6 detectors. The top graph 
gives the angular error as a function of a and sin((^).For clarity, the values above 1° have been set 
to l°.The bottom plots presents the angular error distribution and the fraction of events below a 
given angular threshold. 

leading to a false alarm rate about 10^^ Hz when performing a threefold coincidence). This 
condition is satisfied for 60 % of the sky and the median angular error in this case is 1.3°. As 
a resolution of 1° is obtained for 30 % of the events satisfying the SNR condition, it means 
that about 20 % of the whole sky is seen with an angular error lower than 1°. 

Adding other gravitational waves detectors allows to reduce the blind regions and to lower 
the mean resolution. In the best considered case (6 detectors), the resolution is about 0.4° 
and 99% of the sky is seen with a resolution lower than 0.8°. 

All quoted resolutions (about one degree) are similar to those dehvered by 7-ray satellites 
when the first GRB counterparts have been identified. So, we can expect it will be also 
sufficient for the first identification of gravitational wave sources. 
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